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Abstract 

Simulation calculations are conducted for Si substrates to analyze formation 
and diffusion energies of interstitials under strain condition using statics methods 
based on a Stillinger- Weber type potential function. Defects in the vicinity of the 
surface region and in the bulk are examined, and the role played by compressive and 
tensile strains on the energetics of interstitials is investigated. Results indicate that 
strain alters defect energetics which, in turn, modifies their diffusion characteristics. 
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Introduction 


Diffusion of point defects plays an important role in semiconductor device fab- 
rication technology, so that an understanding of diffusion characteristics of inter- 
stitials in semiconducting materials under strain is highly desireable [1], Presently, 
little is known about the effect of strain on point defect migration properties [2], The 
ability to control strain states may provide additional possibilities for controlling 
diffusion related phenomena in semiconductors. Recently, in several experimental 
studies the role played by strained surface layers on nucleation processes of defects 
has been investigated [3]. However, there is still no consensus about the real mech- 
anism for the formation of defects within the surface region and the effect of strain 
upon it. Experiments on intrinsic point defects have produced conflicting results 

[4]- 


In this work we investigated the effect of strain on the energetics and diffusion 
characteristics of Si interstitials. Furthermore, calculations are conducted to analyze 
variations in the interstitial formation energy and in the diffusion barrier height as 
a function of distance from the exposed surface. 

Calculations 

Energy calculations were carried out employing the empirical potential function 
developed by Stillinger and Weber for systems containing Si species. This function 
has been used successfully in various occasions producing acceptable results [5,6]. 
Calculations were performed employing a molecular statics method representing the 
low temperature limit, T=^0 K. The computational cell used in this study consisted 
of 16 layers containing 16 Si atoms per layer. The bottom two layers were held 
rigid and a two-dimensional periodic boundary condition, PBC, was imposed for 
continuity in simulations involving surface regions. An exposed Si(100)-(2xl) sur- 
face was considered in all computations. A three-dimensional PBC was applied for 
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simulating bulk configurations. In the computational cell the z-direction is taken 
perpendicular to the exposed (100)-(2xl) surface, while the x- and y-directions are 
parallel and perpendicular to the surface dimer rows, respectively. In all cases, the 
diffusion of the interstitial atom is taken parallel to the x-direction (this represents 
the 11 direction along the dimer rows). Strain conditions were imposed unilater- 
ally by elongating or compressing the bulk lattice in the x-, y- and z-directions. 
Calculations include up to 2% compression and up to 3% tension. 

Interstitial formation energies are calculated as 

P n + 1 F 

tf = ■C'n+l ' "n 

n 

where, e f denotes the interstitial formation energy, E n+i represents the total op- 
timized energy of the system containing n + 1 particles, including the interstitial 
atom, and E n is the total optimized energy of defect free system with n atoms. The 
barrier height, A E, along the diffusion path is calculated as [7] 

AE = e f (S)-e f (T) 

Here, e/(S) and e/(T) represent interstitial energies at the saddle point (S) and 
at the tetragonal site (T), respectively. The Stillinger- Weber potential used in 
this study produces the tetragonal site as the energetically more favorable in- 
terstitial site in the Si crystal [5]. The hexagonal site (H), in this case, co- 
incides with the saddle point along the diffusion path, which takes place via 
tetragonal— ►hexagonal— ►tetragonal sites in the crystal. Figure 1 depicts 3D rep- 
resentations of tetragonal and hexagonal sites. 

Per cent variations in energy values, Pe , due to applied strains are calculated 
as: 

P E = 100 • (e,/ej - 1) 

where, e f and ej represent interstitial energies for strained and unstrained systems, 
respectively. In a si mil ar way, the per cent variation in the barrier height, Pe(B) 
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is estimated as: 


Pe{B) = 100 • ( AE/AE° - 1) 

with A E and A E° denoting barrier heights for strained and unstrained cases. 

Figure 2 shows a side view ( xz plane) of the Si crystal depicting the T and H 
sites, and the diffusion path in the general x direction. This path represents the 
energetically lowest path and it was determined by moving the interstitial atom in 
the x direction with incremental steps from its equilibrium (T) site. After every 
displacement, the system was completely equilibrated provided that the moving 
interstitial atom was forced to remain within the plane perpendicular to the diffusion 
path. The path is a straight line for the T-H-T portion which makes an angle of 
about 35 degrees angle with the x-axis. All the T and H sites depicted in Figure 2 
are in the same xz plane for the diffusing interstital along the x-direction. 

Results and Discussions 

Figure 3 shows the per cent variation in the formation energy of a tetragonal- 
site interstitial as a function of applied strain. These results indicate that the 
interstitial formation process is more favorable under tension and it becomes ener- 
getically less favorable under compression. The percentage variation in the intersti- 
tal energy is alm ost linear for strains applied in the z-direction. Due to sy mm etry, 
both tension and compression applied in x- and y-directions show an almost iden- 
tical variation in Pe. The general trend that we obtained here is consistent with 
Antonelli and Bernholc’s results for pressure effects on the formation energy of 
T-type interstitials in Si [8]. 

Variation in the formation energy of an H-site interstitial is shown in Figure 
4 as a function of the applied strain. Similar to the tetragonal site, the interstitial 
formation process is more favorable under tension and it becomes energetically less 
favorable with increasing compression. In this case, as anticipated, variations in the 
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Pe value differ depending on the applied strain directions (i.e. x- or y- directions). 

Per cent variations in calculated energy barrier for interstitial diffusion taking 
place under varying external strains are shown in Figure 5. These results indicate 
that strain applied in the y-direction has almost no effect on AE. In both x- and 
y-directions the applied tension increases Pe(B), impeding the diffusion, and the 
compression lowers Pe(B ), facilitating the interstitial diffusion process. These re- 
sults are consistent with the pressure effect investigations of Antonelli and Bernholc 
[8] on diffusion energies of interstitials in silicon. This outcome is also in agreement 
with the experimental findings of Kringhoj et al [9] on the diffusion of Sb intersti- 
tials, who concluded that diffusivity is enhanced by compressive strain and retarded 
by tensile strain (also see Ref. [1]). Furthermore, experimental findings by Zhao, et 
al. [10] reveal that boron diffusion in silicon is enhanced with increasing pressure; 
these findings also are qualitatively in agreement with present results. 

Relative interstitial formation energies (for T-sites) are depicted in Figure 6 as 
a function of the distance, d z from the exposed surface. Results indicate that e/ 
approaches its bulk value for interstitials deeper than d 2 « 10 A. For d z values be- 
tween approximately 5 to 10 A, the deviation in ef(z)/ef(bulk) is rather small. For 
d z values less than 5 A, however, the relative interstitial formation energy dimin- 
ishes considerably as the interstitial approaches to the exposed surface region. This 
is partially because of the surface reconstruction which alters the overall geometry 
within the top surface layers. Present results, in agreement with earlier work [5], 
indicate that in addition to the dimerization process of top layer atoms, consider- 
able layer-by-layer displacements take place within the top several surface layers. 
The first interlayer spacing (which is the distance between the top and the second 
layer atoms) contracts more than 8% with respect to bulk interlayer spacing. The 
second interlayer spacing, on the other hand, expands about 1.4%, and the third 
interlayer spacing contracts somewhat less than 1%. As a result, within the first 
and second interlayer spacings in particular, geometric features of T and H sites 
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are distorted. Therefore, the interstitial formation energies for these sites deviate 
considerably from their bulk values, which in turn affect diffusion characteristics of 
interstitials. 

Variations in the barrier height (with respect to the bulk value) are shown in 
Figure 7 as a function of d z . Calculated A E values for diffusion display a steady 
increase as the defect approaching the exposed surface region. This outcome indi- 
cates that interstitial diffusion within the top few layers is not favored energetically. 
Values for the top three interlayer spacings are not shown in the graph. It appears 
that the diffusion process via T-*H-+T sites within these very top layers is not 
energetically feasible. Probably, interstitial motion in these cases takes place via 
different diffusion mechanisms which we are presently investigating. 

Conclusions 

Formation energies of interstitials and their diffusion characteristics were ana- 
lyzed for silicon systems under uniaxial strain conditions. For both T and H sites, 
the interstitial formation process is energetically more favorable under tension but 
less favorable under compressive strain. The trend is s imil ar for strains in x-, y- 
and z-directions. Results indicate that strain applied in the direction perpendicu- 
lar to a diffusing interstitial (y-direction) has a minimal or no effect on A E, the 
barrier height. In x- and z-directions, on the other hand, tension increases AE, 
and compression was found to reduce A E and thus to facilite the diffusion process. 
Furthermore, present calculations show that the relative formation energy of an 
interstitial decreases within the surface region. This is partially because of the con- 
figurational changes (such as dimerization and the multilayer relaxation) which take 
place within the exposed surface region. At the same time, however, the barrier 
height for the interstital diffusion process via tetragonal— ►hexagonal— ►tetragonal 
sites increases as the defect nears the exposed surface region. 
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Figure Captions 

Figure 1. Three-dimensional views of (a) tetragonal and (b) hexagonal sites. The 
darker sphere represents the interstitial Si atom, and bonded light color spheres are 
the lattice Si atoms. 

Figure 2. Schematic side view representation of the Si crystal depicting the xz 
plane. The dotted line shows the diffusion path. Open and solid circles represent 
tetrahedral and hexagonal sites, respectively. 

Figure 3. Per cent variation in energy, P E , f° r a tetragonal-site interstitial versus 
the applied strain. Solid, dashed and dotted lines show variations in P E for strain 
applied in the x-, y- and z-directions, respectively. Negative and positive values of 
strain represent compression and tension, respectively. 

Figure 4. Per cent variation in energy, P E , for a hexagonal-site interstitial versus 
the applied strain. Solid, dashed and dotted lines show variations in P E for strain 
applied in the x-, y- and z-directions, respectively. Negative and positive values of 
strain represent compression and tension, respectively. 

Figure 5. Per cent variation in the barrier height, P E (B), for an interstitial diffu- 
sion as a function of applied strain. Solid, dashed and dotted lines show variations 
in P E {B) for strain applied in the x-, y- and z-directions, respectively. Negative 
and positive values of strain represent compression and tension, respectively. 

Figure 0. Variations in the interstitial formation energy as a function of the 
distance, d z , from the exposed surface. Relative energy values with respect to bulk 
are plotted on the ordinate versus d z . The values of d z are given in A. 

Figure 7. Variations in the diffusion barrier height as a function of the distance, 
d z , from the exposed surface. Relative values of the barrier heights with respect to 
the bulk value are plotted versus d z . The values of d z are given in A. 
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